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Abstract 

We highlight issues which are often underestimated in the experimental analyses on 
quarkonium polarization: the relation between the parameters of the angular distri- 
butions and the angular momentum composition of the quarkonium, the importance 
of the choice of the reference frame, the interplay between observed decay and produc- 
tion kinematics, and the consequent influence of the experimental acceptance on the 
comparison between experimental measurements and theoretical calculations. Given 
the puzzles raised by the available experimental results, new measurements must pro- 
vide more detailed information, such that physical conclusions can be derived without 
relying on model-dependent assumptions. We describe a frame-invariant formalism 
which minimizes the dependence of the measurements on the experimental accep- 
tance, facilitates the comparison with theoretical calculations, and probes systematic 
effects due to experimental biases. This formalism is a direct and generic consequence 
of the rotational invariance of the dilepton decay distribution and is independent of 
any assumptions speciflc to particular models of quarkonium production. The use 
of this improved approach, which exploits the intrinsic multidimensionality of the 
problem, will significantly contribute to a faster progress in our understanding of 
quarkonium production, especially if adopted as a common analysis framework by 
the LHC experiments, which will soon perform analyses of quarkonium polarization 
in proton-proton collisions. 
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1 Introduction 



Detailed studies of quarkonium prodution should provide significant progress in our 
understanding of quantum chromodynamics (QCD) [IJ. However, our present under- 
standing of this physics topic is rather limited, despite the multitude of experimental 
data accumulated over more than 30 years. The differential ^/ip and %p' direct 
production cross sections measured (in the mid 1990's) by CDF, in pp collisions at 
1.8 TeV [2], were seen to be around 50 times larger than the available expectations, 
based on leading order calculations made in the scope of the Colour Singlet Model. 
The non-relativistic QCD (NRQCD) framework [3], where quarkonia can also be pro- 
duced as coloured quark pairs, succeeded in describing the measurements, opening a 
new chapter in the studies of quarkonium production physics. However, these calcu- 
lations depend on non-perturbative parameters, the long distance colour octet matrix 
elements, which have been freely adjusted to the data, thereby decreasing the impact 
of the resulting agreement between data and calculations. More recently, calcula- 
tions of next-to-leading-order (NLO) QCD corrections to colour-singlet quarkonium 
production showed an important increase of the high-px rate, significantly decreasing 
the colour-octet component needed to reproduce the quarkonium production cross 
sections measured at the Tevatron PJ. 

Given this situation, differential cross sections are clearly insufficient information 
to ensure further progress in our understanding of quarkonium production. Experi- 
mental studies of the polarization of the J^^ = 1 quarkonium states, which decay 
into lepton pairs, will certainly provide very useful complementary information. In 
fact, the competing mechanisms dominating in the different theoretical approaches 
lead to very different expected polarizations of the produced quarkonia. On one hand, 
the NRQCD calculations [SIEIIZI, dominated by the colour-octet component, predict 
that, at Tevatron or LHC energies and at asymptotically high px, the directly pro- 
duced ip' and ^/ip mesons are produced almost fully transversely polarized (i.e. with 
dominant angular momentum component Jz = ±1) with respect to their own mo- 
mentum direction (the helicity frame) . On the other hand, according to the new NLO 
calculations of colour-singlet quarkonium production ^ these states should show a 
strong longitudinal (J^ = 0) polarization component. 

Having two very different theoretical predictions appears to be an ideal situation 
when seen from an experimentalist's perspective, as one may think that it should be 
relatively straightforward to discriminate between the two theory frameworks using 
experimental measurements. Somewhat surprisingly, however, this is not the case. In 
fact, the present experimental knowledge is incomplete and contradictory. Studies of 
the ip' polarization have been published on the basis of data collected by the CDF II 
experiment [8]. Unfortunately, the large experimental uncertainties caused by the 
small size of the data samples prevent from drawing meaningful conclusions. In 
principle, a more precise test of the theoretical predictions should be provided by the 
J/ip data, given their much higher statistical accuracy. However, the experimental 
perspective is more complicated in this case, because a significant fraction (around 
one third [9]) of promptly produced J/tp mesons (i.e. excluding contributions from 
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B hadron decays) comes from Xc and ip' feed-down decays. This sizeable source of 
indirectly produced i/ip mesons is not subtracted from the current measurements, 
and its kinematic dependence is not precisely known. Despite this limitation, it seems 
safe to say that the pattern measured by CDF [8J of a slightly longitudinal polarization 
of the inclusive prompt J/^'s is incompatible with any of the two theory approaches 
mentioned above. The situation is further complicated by the intriguing lack of 
continuity between fixed-target and collider results, which can only be interpreted in 
the framework of some specific (and speculative) assumptions still to be tested |10j . 

The hh system should satisfy the non-relativistic approximation much better than 
the cc case. For this reason, the T data are expected to represent the most decisive test 
of NRQCD. However, the comparison with the existing T(IS') polarization data from 
Tevatron [II1|12] is far from conclusive. The results indicate that, for px < 15 GeV/c, 
the T(IS') is produced either unpolarized (CDF) or longitudinally polarized (DO) in 
the helicity frame, and this discrepancy cannot be reasonably attributed to the differ- 
ent rapidity windows covered by the two experiments. Furthermore, the precision of 
the data for pi higher than 15 GeV/c is not sufficient to provide a significant test of 
the crucial hypothesis that very high-px quarkonia, produced by the fragmentation of 
an outgoing (almost on-shell) gluon, are fully transversely polarized along their own 
direction. At lower energy and pt, the E866 experiment [13] has shown yet a different 
polarization pattern: the T(25') and T(3S') states have maximal transverse polariza- 
tion, with no significant dependence on transverse or longitudinal momentum, with 
respect to the direction of motion of the colliding hadrons (Collins-Soper frame). Un- 
expectedly, the T(IS'), whose spin and angular momentum properties are identical 
to the ones of the heavier T states, is, instead, found to be only weakly polarized. 
These results give interesting physical indications. First, the maximal polarization of 
T(2S') and T(3S') along the direction of the interacting particles places strong con- 
straints on the topology and spin properties of the underlying elementary production 
process. Second, the small T(IS') polarization suggests that the bottomonium family 
may have a peculiar feed-down hierarchy, with a very significant fraction of the lower 
mass state being produced indirectly; at the same time, the polarization of the T's 
coming from Xb decays should be substantially different from the polarization of the 
directly produced ones. 

This rather confusing situation demands a significant improvement in the accu- 
racy and detail of the polarization measurements, ideally distinguishing between the 
properties of directly and indirectly produced states. We remind that the lack of a 
consistent description of the polarization properties represents today's biggest uncer- 
tainty in the simulation of the LHC quarkonium production measurements and will 
probably be the largest contribution to the systematic error affecting the measure- 
ments of quarkonium production cross sections and kinematic distributions. Indeed, 
the probability that the detector accepts lepton pairs resulting from decays of quarko- 
nium states is strongly dependent on the polarization of those states. Therefore, even 
from a purely experimental point of view it is very unsatisfactory that essential prop- 
erties of these objects, such as kinematic details of how they decay into lepton pairs 
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(on which their reconstruction is based), are subjected to such a high degree of un- 
certainty. 

It is true that measurements of the quarkonium decay angular distributions are 
challenging, multi-dimensional kinematic problems, which require large event samples 
and a very high level of accuracy in the subtraction of spurious kinematic correlations 
induced by the detector acceptance. The complexity of the experimental problems 
which have to be faced in the polarization measurements is testified, for example, 
by the disagreement between the CDF results obtained in Run I and Run 11 for 
the ^/ip [m [8] and by the contradictory results obtained by CDF and DO for the 
T(15'). However, it is also true, as we shall emphasize in this paper, that most exper- 
iments have exploited, and presented in the published reports, only a fraction of the 
physical information derivable from the data. This happens, for example, when the 
measurement is performed in only one polarization frame and is limited to the polar 
projection of the decay angular distribution. As we have already argued in Ref. [TO] , 
these incomplete measurements do not allow definite physical conclusions. At best, 
they confine such conclusions to a genuinely model-dependent framework. Moreover, 
such a fragmentary description of the observed physical process obviously reduces the 
chances of detecting possible biases induced by not fully controlled systematic effects. 

In this paper we review the mathematical framework for the description of the 
observable polarization of quarkonium states decaying electromagnetically into lepton 
pairs. We focus our attention on aspects that need to be taken in consideration in the 
analysis of the data, so as to maximize the physical significance of the measurement 
and provide all elements for its unambiguous interpretation within any theoretical 
framework. By increasing the level of detail of the physical information extracted 
from the data, the proposed methodologies also offer the possibility of performing 
consistency checks which can expose unaccounted detector or analysis biases. The 
only relevant theoretical ingredients of our discussion are the quantum properties of 
angular momentum and basic conservation rules of the electromagnetic interaction 
(parity, fermion chirality). All the results presented here are, therefore, valid in 
general for any quarkonium production mechanism. 

In Section |2] we define the concept of polarization and give simple examples of 
how basic production mechanisms can lead to the formation of polarized quarkonium 
states. We then focus on the dilepton decay distribution of "^Si quarkonia, a rela- 
tively simple case, and provide detailed geometric and kinematic considerations. In 
Section |3] we recall the basic principles leading to the general expression describing 
the angular distribution of the decay products, while in Section |4] we describe how 
the observed anisotropy parameters depend on the choice of the reference frame. 
Section |5] is devoted to a detailed description of how the production kinematics in- 
fluences the observed polarization, depending on the quarkonium momentum and on 
the observation frame. We also discuss quantitatively the influence of the intrinsic 
parton transverse momenta on the polarization measurement when the natural axis 
is along the relative flight direction of the colliding partons. In Section [6] we derive 
the existence of a frame-independent identity which relates the observable parameters 
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of the decay distribution to one frame-invariant polarization parameter. We discuss 
how this relation, formally including the Lam- Tung identity [15] as a particular case, 
improves the representation of polarization results and can be used to perform consis- 
tency checks in the experimental analyses. We continue with some remarks, given in 
Section [7], on how the existence of intrinsic parton transverse momentum affects the 
polarization measurement. We conclude, in Section [8} with a few examples inspired 
from existing experimental measurements, which should provide concrete evidence 
for the usefulness of the approaches discussed in this paper, in view of ensuring an 
improved understanding of quarkonium production. 

2 Basic polarization concepts 

Because of angular momentum conservation and basic symmetries of the electro- 
magnetic and strong interactions, a particle produced in a certain superposition of 
elementary mechanisms may be observed preferentially in a state belonging to a defi- 
nite subset of the possible eigenstates of the angular momentum component along 
a characteristic quantization axis. When this happens, the particle is said to be po- 
larized. Figure [T] shows examples of leading-order diagrams of elementary production 
processes giving rise to different types of polarizations. 




Figure 1: Examples of leading-order diagrams for production mechanisms giving rise 
to observable polarizations: (a) vector quarkonium production in electron-positron 
annihilation; (b) Drell-Yan production in quark-antiquark annihilation; (c) quarko- 
nium production by gluon fragmentation to colour-octet cc. 

Vector (J = 1 ) quarkonia have the same charge-parity as an electron-positron 
pair and can be produced in electron-positron annihilation, via an intermediate pho- 
ton (Fig.[l]a). The states originating from this process are polarized, as a consequence 
of helicity conservation, a general property of QED in the relativistic (massless) limit. 
The dynamics of the coupling of electrons to photons is, in fact, described by terms 
of the form u'-y^u = Ui^'-f'^ui^ + mr7'^mr, where are the Dirac matrices, u is the 
electron spinor, and L (R) indicate its left-handed (right-handed) chiral components. 
Terms with opposite chiral components are absent, meaning that the fermion chiral- 
ity is preserved in the interaction with a photon. When the fermions are assumed 
to have zero mass, so that the direction of their momenta cannot be reversed by 
any Lorentz transformation, left-handed and right-handed chiral components become 
eigenstates of the helicity operator h = S ■ p/\p\, corresponding to the projection of 
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the spin on the momentum direction. In this case, chirahty conservation becomes he- 
hcity conservation. In the diagram of Fig. [T]a, this rule implies that the annihilating 
electron and positron must have opposite helicities, because the intermediate photon 
has zero (fermion) helicity. Since in the laboratory their momenta are opposite, their 
spins must be parallel. Because of angular momentum conservation, the produced 
quarkonium has, thus, angular momentum component = ±1 along the direction 
of the colliding leptons. This precise QED prediction (the relative amplitude for the 
Jz = component is of order rrie/Ee ^ 3 x 10^^ for J ftp production and smaller for T 
production) is commonly used as a base assumption in quarkonium measurements in 
electron-positron annihilations (as, for example, in the recent analysis of Ref. [IB])- 
The fact that the dilepton system coupled to a photon is a pure = ±1 state is also 
an essential ingredient in the determination of the expression for the dilepton-decay 
angular distributions of vector quarkonia (see Section |3]). 

The same reasoning can be applied to the production of Drell-Yan lepton pairs 
in quark-antiquark annihilation (Fig. [l]b): the quark and antiquark, in the limit of 
vanishing masses, must annihilate with opposite helicities, resulting in a dilepton 
state having = ±1 along the direction of their relative velocity. The experimental 
verification of this basic mechanism has reached an impressive level of accuracy [13] . 
Quark helicity is conserved also in QCD, when the masses can be neglected. Similarly 
to the Drell-Yan case, quarkonia originating from quark-antiquark annihilation (into 
intermediate gluons) will thus tend, provided they are produced alone, to have their 
angular momentum vectors "aligned" (J^ = ±1) along the beam direction. This 
prediction is in good agreement with the xa, Xc2 and ip' polarizations measured in 
low-energy proton-antiproton collisions [TTj ITS]. 

At very high px, quarkonium production at hadron colliders should mainly pro- 
ceed by gluon fragmentation (TH]. In NRQCD, heavy-quark velocity scaling rules 
for the non perturbative matrix elements, combined with the as and 1/pt power 
counting rules for the parton cross sections, predict that J/ip and ip' production 
at high pt is dominated by gluon fragmentation into the color-octet state ccpS"!^^] 
(Fig. [l]c). Transitions of the gluon to other allowed colour and angular momentum 
configurations, containing the cc in either a colour-singlet or a colour-octet state, 
with spin S* = 0, 1 and angular momentum L = 0,1,2,..., as well as additional 
gluons (cc[^5'q^^]c/, cc['^Pf'^]g, cc['^S[^^]gg, etc.), are more and more suppressed with 
increasing p^. Up to small corrections, the fragmenting gluon is believed to be on 
shell and have, therefore, helicity ±1. This property is inherited by the ccpsj^^] 
state and remains intact during the non-perturbative transition to the colour-neutral 
physical state, via soft-gluon emission. In this model, the observed charmonium has, 
thus, angular momentum component Jz = ±1, this time not along the direction of 
the beam, but along its own flight direction. 

"Unpolarized" quarkonium has the same probability, 1/(2 J + 1), to be found in 
each of the angular momentum eigenstates, Jz = —J, — J + 1, . . . , + J. This is the 
case, for example, in the colour evaporation model [20]. In this framework, similarly 
to NRQCD, the QQ pair is produced at short distances in any colour and angular 
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momentum configuration. However, contrary to NRQCD, no hierarchy constraints 
are imposed on these configurations, so that the cross section turns out to be domi- 
nated by QQ pairs with vanishing angular momentum (^5'o), in either colour-singlet 
or colour-octet states. In their long distance evolution through soft gluon emissions, 
J = states get their colour randomized, assuming the correct quantum numbers of 
the physical quarkonium. As a result, the final angular momentum vector J has no 
preferred alignment. 

In two-body decays (such as the '^5'i — )■ case considered in this paper), the 
geometrical shape of the angular distribution of the two decay products (emitted back- 
to-back in the quarkonium rest frame) refiects the polarization of the quarkonium 
state. A spherically symmetric distribution would mean that the quarkonium would 
be, on average, unpolarized. Anisotropic distributions signal polarized production. 




Figure 2: The coordinate system for the measurement of a two-body decay angular 
distribution in the quarkonium rest frame. The y axis is perpendicular to the plane 
containing the momenta of the colliding beams. The polarization axis z is chosen 
according to one of the possible conventions described in Fig. |3j 

The measurement of the distribution requires the choice of a coordinate system, 
with respect to which the momentum of one of the two decay products is expressed 
in spherical coordinates. In inclusive quarkonium measurements, the axes of the 
coordinate system are fixed with respect to the physical reference provided by the 
directions of the two colliding beams as seen from the quarkonium rest frame. Figure[2] 
illustrates the definitions of the polar angle determined by the direction of one of the 
two decay products (e.g. the positive lepton) with respect to the chosen polar axis, and 
of the azimuthal angle (p, measured with respect to the plane containing the momenta 
of the colliding beams ("production plane"). The actual definition of the decay 
reference frame with respect to the beam directions is not unique. Measurements 
of the quarkonium decay distributions have used three different conventions for the 
orientation of the polar axis (see Fig. [3]): the direction of the momentum of one of the 
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production plane y 




Figure 3: Illustration of three different definitions of the polarization axis z (CS: 
Collins-Soper, GJ: Gottfried- Jackson, HX: helicity) with respect to the directions of 
motion of the colliding beams (bi, b2) and of the quarkonium (Q). 

two colliding beams (Gottfried- Jackson frame |2T], GJ), the opposite of the direction 
of motion of the interaction point (i.e. the flight direction of the quarkonium itself 
in the center-of-mass of the colliding beams: helicity frame, HX) and the bisector 
of the angle between one beam and the opposite of the other beam (Collins-Soper 
frame [22j, CS). The motivation of this latter definition is that, in hadronic collisions, 
it coincides with the direction of the relative motion of the colliding partons, when 
their transverse momenta are neglected (the validity and limits of this approximation 
are discussed in detail in Section [?]). For our considerations, we will take the HX 
and CS frames as two extreme (physically relevant) cases, given that the GJ polar 
axis represents an intermediate situation. We note that these two frames differ by a 
rotation of 90° around the y axis when the quarkonium is produced at high and 
negligible longitudinal momentum (px ^ |pl|)- All definitions become coincident in 
the limit of zero quarkonium p^. In this limit, moreover, for symmetry reasons any 
azimuthal dependence of the decay distribution is physically forbidden. 

We conclude this section by defining the somewhat misleading nomenclature which 
is commonly used (and adopted, for convenience, also in this paper) for the polar- 
ization of vector mesons. These particles share the quantum numbers of the photon 
and are therefore said, by analogy with the photon, to be "transversely" polarized 
when they have spin projection = ±1. The counterintuitive adjective originally 
refers to the fact that the electromagnetic field carried by the photon oscillates in 
the transverse plane with respect to the photon momentum, while the photon spin is 
aligned along the momentum. "Longitudinal" polarization means Jz = 0. By further 
extension, the same terms are also used to describe the "spin alignment" of vector 
quarkonia not only with respect to their own momenta (HX frame), but also with 
respect to any other chosen reference direction (such as the GJ or CS axes). 
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3 Dilepton decay angular distribution 



Vector quarkonia, such as the J/?/;, ip' and T{nS) states, can decay electromagneti- 
cally into two leptons. The reconstruction of this channel represents the cleanest way, 
both from the experimental and theoretical perspectives, of measuring their produc- 
tion yields and polarizations. In this and the following sections we discuss how to 
determine experimentally the "spin alignment" of a vector quarkonium by measuring 
the dilepton decay angular distribution. For convenience we mention explicitly the 
J/V' as the decaying particle, but considerations and results are valid for any J = 1 
state. 

We begin by studying the case of a single production "subprocess" , here defined 
as a process where the J/ip is formed as a given superposition of the three J = 1 
eigenstates, Jz = +1, —1, with respect to the polarization axis z: 

|\/) = 6+i| + l)+6_i|-l)+6o|0). (1) 

The calculations are performed in the J/ip rest frame, where the common direction of 
the two leptons define the reference axis z', oriented conventionally along the direction 
of the positive lepton. The adopted notations for axes, angles and angular momentum 
states are illustrated in Fig. [4j Because of helicity conservation for (massless) fermions 




Figure 4: Sketch of the decay i /ip ^ , showing the notations used in the text 
for axes, angles and angular momentum states. 

in QED, the dilepton system coupled to a photon in the process J/V' — 7* — has 
angular momentum projection ±1 along z', i.e. it can be represented as an eigenstate 
of Jz', \(-^(-~] 1, V) with /' = +1 or — 1. We want to express this state as a superposition 
of eigenstates of J^, 1, /) with / = 0, ±1. To perform this change of quantization 

axis, we use a general result of angular momentum theory, which we recall in the 
following paragraphs. 

We indicate by /3, 7) the rotation from a generic set of axes (x, z) to the set 
{x',y', z'), a, (3 and 7 denoting the Euler angles. Positive rotations are defined by the 
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right-hand rule. An eigenstate | J, M') of J^/ can then be expressed as a superposition 
of the eigenstates | J, M) of Jz through the rotation transformation 



+J 



|J,M')= ^ Vi,,{R)\J,M). (2) 

M=-J 

The (complex) rotation matrix elements 'D'ljj^, are defined as 

HiM'icy, (3, 7) = e-^''-dUAP)e-^'''' (3) 
in terms of the (real) reduced matrix elements 

Tam(J+M,J-M') 
t=max(0,M-M') 

v/( J + M)! ( J - M)! ( J + MQ! (J - MQ! 
^ ( J + M - t)! ( J - M' - t)! t! (t - M + M')! ^ 

p\2J+M-M'-2t / p\ 2t-M+M' 



xicos-j (^sm^ 

The rotation we need in our case has the effect of bringing one quantization axis (z) 
to coincide with another {z'). The most general rotation performing this projection 
can be parametrized with P = d and a = —7 = ip. The dilepton angular momentum 
state is therefore expressed in terms of eigenstates of Jz as 

\tr;i,i')= E vUif,^,-^)\e+e-;i,i). (5) 

l=0,±l 



The amplitude of the partial process J/tp{m) — )• i'^i (/') represented in Fig. |4]is 



«=0,±1 
\1* 



= BVl:,{^,^,-if), (6) 

where we imposed that the transition operator B is of the form 
\B\ J/ip]l,m) = B 5ml because of angular momentum conservation, 
with B independent of m (for rotational invariance). The total amplitude for 
J/ip ^ where the J/ip is given by the superposition written in Eq. [l| is 



B, = bmBVl:i,i^,{},-^) 

m=-l,+l 

= J2 CirnV]:,{^,d,-^). (7) 



m=-l,+l 
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The probability of the transition is obtained by squaring Eq. [7]and summing over the 
(unobserved) spin ahgnments (/' = ±1) of the dilepton system, with equal weights 
attributed, for parity conservation, to the two configurations. Using Eq. [3| with 
c?o ±1 = ± sin ^9/^2, = (1 + cos'^)/2 and = (1 - cosi9)/2 (from Eq. |3|, 

one obtains the angular distribution 

W{cos'd,ip) oc \Bi>\'^ oc (1 + X^cos^'d 

itii ^3 + 

+ sin^ cos 2ip + A^<^ sin 2{} cos ip (8) 

+ sin^ sin 2(p + A;^^ sin 2'i9 sin ip) , 

with M = |aop + + and 

Ai9 



Xip 



A/'-3 




2 






2 



A 



2Re[a5i*a_i] 

y2Re[4'^*(a+i -a_i)] 

A/'+ kol 
— 2Im[a!i.ia_i] 



2 



(9) 



A 



^ ^f+\ao\^ ' 

_ -A/2Im[ao(a+i + a_i)] 



Figure [5] shows the shapes of the dilepton decay distributions in the two polarization 
cases m = ±1 (a) and m = (b); the m = +1 and m = —1 configurations are indis- 
tinguishable because of rotation invariance. The same distributions are also shown, 
in panels (c) and (d), as seen when studied in frames rotated by 90°, anticipating the 
discussion in Section |4l 

It is worth noticing that it is impossible to chose the decay amplitudes am and, 
therefore, the component amplitudes bm such that all decay parameters in Eq. [8] 
vanish. This means that the angular distribution of the decay of a J = 1 state is 
never intrinsically isotropic. Even if it is conceivable that a lucky superposition of 
different production processes might lead to a fortuitous cancellation of all decay 
parameters, such an exceptional case would signal a non-trivial physical polarization 
scenario, caused by spin randomization effects, or (semi-)exclusive configurations in 
which the observed state is produced together with certain final state objects. In other 
words, polarization is an essential property of the quarkonium states. This remark is 
particularly relevant when we consider that all existing Monte Carlo generators use 
an isotropic dilepton distribution as the default option for quarkonium production in 
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a) 



X 



m = ±l 



b) 



m-0 



c) 



y' 




Figure 5: Graphical representation of the dilepton decay distribution of "trans- 
versely" (a) and "longitudinally" (b) polarized quarkonium in the natural frame and 
in frames rotated by 90° (c-d). The probability of the lepton emission in one direction 
is represented by the distance of the corresponding surface point from the origin. 



hadronic collisions, a non-trivial assumption with a strong influence on the acceptance 
estimates and, therefore, on both normalizations and kinematic dependencies of the 
measured quarkonium cross sections. 

In this paper we only consider inclusive production. Therefore, the only possible 
experimental definition of the xz plane coincides with the production plane, contain- 
ing the directions of the colliding particles and of the decaying particle itself. The 
last two terms in Eq. [8] introduce an asymmetry of the distribution by reflection 
with respect to the production plane, an asymmetry which is not forbidden in indi- 
vidual (parity-conserving) events. In hadronic collisions, due to the intrinsic parton 
transverse momenta, for example, the "natural" polarization plane does not coincide 
event-by-event with the experimental production plane. However, the symmetry by 
reflection must be a property of the observed event distribution when only parity- 
conserving processes contribute. Indeed, the terms in sin^ sin 2ip and sin 2^9 sin ip are 
unobservable, because they vanish on average. 

In the presence of n contributing production processes with weights /'•*•', the most 
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general observable distribution can be written as 



+ \ siv? 'd cos 2ip + A^^ sin 2'd cos ip) , 

where W^"^^ (cos "(9, is the "elementary" decay distribution corresponding to a single 
subprocess (given by Eqs. |8]and|9| adding the index [i] to the decay parameters) and 
each of the three observable shape parameters, X = A^, A^ and A^^, is a weighted 
average of the corresponding parameters, X^^\ characterizing the single subprocesses, 

trS + AW /^3 + A« 

We conclude this section with the derivation of formulae which can be used for 
the determination of the parameters of the observed angular distribution, as an alter- 
native to a multi-parameter fit to the function in Eq. 10 The integration over either 
or cost? leads to one-dimensional angular distributions, 

W{cos^) oc (l + A^cos^t?) , (12) 

3 + A^ 

2A 

W{v) oc l + --^cos2^, (13) 

from which A^ and A^ can be determined in two separate steps, possibly improving 
the stability of the fit procedures in low-statistics analyses. The "diagonal" term, 
A^(^, vanishes in both integrations but can be extracted, for example, by defining the 
variable (f as 

r (/.-fTT for cos^<0 
^ \ y^-l for cos^9>0 ^ ^ 

(adding or subtracting 27r when ip does not fall into one continuous range, e.g. [0, 27r]) 
and measuring the distribution 

W{^) oc 1 + ^-^ cos^. (15) 
3 + A^ 

Each of the three parameters can also be expressed in terms of an asymmetry be- 
tween the populations of two angular topologies (which are equiprobable only in the 
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unpolarized case): 



(16) 



P(|cos^9| > 1/2) -P(|cos?9| < 1/2) _ 3 A^^ 
P(| cosd] > 1/2) + P(| cos^9| < 1/2) ~ 4 3 + ' 

P(cos2(^ > 0) - P(cos2v9 < 0) _ 2 2A^ 
P(cos2(^ > 0) + P(cos2(^ < 0) ~ ^ 3 + A^ ' 

P(sin2?9cos^ > 0) - P(sin2i9cos^ < 0) _ 2 2A^^^ 
P(sin2?9cosv? > 0) + P(sin2'i9cosv2 < 0) ~ tt 3 + A^j ' 

In analyses applying efficiency corrections to the reconstructed angular spectra, the 
use of these formulae may require an iterative re-weighting of the Monte Carlo data, 
in order to compensate for the effect of the non-uniformity of those experimental 
corrections. In "ideal" experiments with uniform acceptance and efficiencies over 
cos '(9 and if (such as in Monte Carlo studies at the generation level) the parameters 
can be obtained from the average values of certain angular distributions: 

2 q\ _ 1 + |Ai) 



(cos' ^) 



3 + A^ ' 



(cos 2^) = (17) 



(sin 219 cosip) 



4 A, 



5 3 + A. 



4 Dependence of the measurement on the obser- 
vation frame 

All possible experimentally definable polarization axes in inclusive measurements be- 
long to the production plane (defined in Fig. [3]). We can, therefore, parametrize 
the transformation from an observation frame to another by one angle describing a 
rotation about the y axis. Instead of rotating the angular momentum state vectors, 
we can apply a purely geometrical transformation directly to the observable angular 
distribution. The rotation matrix 

(cos S — sin S\ 
10 (18) 
sin (5 cos (5 J 

brings the old frame to coincide with the new one, the positive sign of S being defined 
by the right-hand rule (we will discuss in Section [s] how the sign of 5 depends in an 
observable way on the conventions chosen for the orientation of the z and y axes, and 
how, specifically, the angle between the HX and CS axes depends on the quarkonium 
production kinematics). The unit vector f = (sin?9cosv5, sin^9sin<^, cos-i?) indicating 
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the lepton direction in the old frame is then expressed as f = Ry^{6)r' as a function 
of the coordinates in the new frame. In particular, 

cos?? = — sin 5 sin COS v?' + cos^cos??'. (19) 



Substituting Eq. [19] into Eq. [10} we obtain the angular distribution in the rotated 
frame: 

W\cos^',^') oc -^(1 +A;cos2^?' (20) 
3 + A,, 



+ A' sin^ i!)' cos 2(f' + A^ sin 2^9' cos if') 



where 



^ A^-3A ^ Ay + A 

i + A ' i + A ' 

Ai)y cos25 - 1 (Aij - Ay) sin25 
= ^ , ^^1) 

with ^ = 2 ^"^^ ~ ^"^^ ^^^"^ ^ ~ 2 ^^"^ ^^'^ ' 

Since the magnitude of the "polar anisotropy", A^, never exceeds 1 in any frame, we 
deduce the frame-independent inequalities 

< ^(1 + A^), lA^yl < ^(1-Ay), (22) 

which imply the bounds |Ay| < 1 and |A^y| < 1. More interestingly, we can see 
that I Ay I < 0.5 when A^ = and must vanish when A^ — —1. The most general 
phase space for the three angular parameters is represented in Fig. [6j There is an 
alternative notation, widespread in the literature, where the coefficients A, v/2 and 
/i replace, respectively, A^, Ay and A^y. In that case, hence, we have |z/| < 2. 

To illustrate the importance of the choice of the observation frame, we consider 
specific examples assuming, for simplicity, that the observation axis is perpendicu- 
lar to the natural axis {5 = ±90°). This case is of physical relevance since when 
the decaying particle is produced with small longitudinal momentum (|pl| ^ Pt, 
a frequent kinematic configuration in collider experiments) the CS and HX frames 
are actually perpendicular to one another. When 6 = 90°, a natural "transverse" 



polarization (A^ = +1 and Ay = A^y = 0), for example, transforms (Eq. 21) into 
an observed polarization of opposite sign (but not fully "longitudinal"), A^ = —1/3, 
with a significant azimuthal anisotropy, A'^ = 1/3, shown in Fig. |5] (c). In terms of 
angular momentum wave functions, a state which is fully "transverse" with respect 
to one quantization axis is a coherent superposition of 50% "transverse" and 50% 
"longitudinal" components with respect to an axis rotated by 90° (Eq. [2]): 

|1,±1) ^ + T ^|1,0). (23) 
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Figure 6: Allowed regions for the decay angular parameters. 



The amplitude of the transition of this mixed state to the "rotated" dilepton state 
in Eq. [s] contains three terms with relative phases (due to the (y? dependence of the 
rotation matrix) giving rise to the observable azimuthal dependence. The same polar 
anisotropy A'^ = —1/3 would be measured in the presence of a mixture of at least 
two different processes resulting in 50% "transverse" and 50% "longitudinal" natural 
polarization along the chosen axis. In this case, however, no azimuthal anisotropy 
would be observed. As a second example, we note that a fully "longitudinal" natural 
polarization (A^ = —1) translates, in a frame rotated by 90° with respect to the 
natural one. Fig. [5](d), into a fully "transverse" polarization (A'^ = +1), accompanied 
by a maximal azimuthal anisotropy (A'^ = — 1). In terms of angular momentum, the 
measurement in the rotated frame is performed on a coherent admixture of states, 

|i,o) ^ 7! '^'+^^"7!'^'-^^' ^^^^ 

while a natural "transverse" polarization would originate from the statistical super- 
position of uncorrelated |1, +1) and |1, —1) states. The two physically very different 
cases of a natural transverse polarization observed in the natural frame, shown in 
Fig. |5] (a), and a natural longitudinal polarization observed in a rotated frame, shown 
in Fig. [5] (d), are experimentally indistinguishable when the azimuthal anisotropy 
parameter is integrated out. These examples show that a measurement (or theoret- 
ical calculation) consisting only in the determination of the polar parameter A^ in 
one frame contains an ambiguity which prevents fundamental (model-independent) 
interpretations of the results. The polarization is only fully determined when both the 
polar and the azimuthal components of the decay distribution are known, or when 
the distribution is analyzed in at least two geometrically complementary frames. 



16 



5 Effect of production kinematics on the observed 
decay kinematics 

Ideally, the dependence of the polarization on the momentum components of the pro- 
duced quarkonium should reflect the relative contribution of individual production 
processes in different kinematic regimes, thereby providing information of fundamen- 
tal physical interest. However, the observations are, in general, affected by some 
experimental limitations, which must be carefully taken in consideration. First, the 
frame-dependent polarization parameters A^, and A^^ can be affected by a strong 



explicit kinematic dependence (encoded in the parameter 6 in Eq. 21), reflecting the 
change in direction of the chosen experimental axis (with respect to the "natural 
axis") as a function of the quarkonium momentum. Second, detector acceptances 
and event samples with limited statistics induce a dependence of the measurement 
on the distribution of events effectively accepted by the experimental apparatus. 

To better explain the first problem, let us consider the HX and CS frames as the 
experimental and natural frames, respectively. We start by calculating the angle be- 
tween the polarization axes of the CS and HX frames as a function of the quarkonium 
momentum. The beam momenta in the "laboratory" frame (centre of mass of the 
colliding particles), written in longitudinal and transverse components with respect 
to the quarkonium direction, are Pi = —P2 = Pcos© ?j| + PsinB i±, where P is 
their modulus and is the angle formed by the quarkonium momentum with respect 
to the beam axis, defined in terms of the quarkonium momentum p as cos = Ph/p- 
When boosted to the quarkonium rest frame, the two vectors become (neglecting 
the masses of the colliding particles) P'l = {jPcosQ — P'jP) i\\ + PsinG i± and 
P'2 = (— 7Pcos0 — P'jP) ill — PsinO where 7 = E/m is the Lorentz factor of 



the quarkonium state, and j3 = p/ E = a/1 — 1/7^- The unit vectors indicating the z 
axis directions in the HX and CS frames are 



P'l + P'2 V 



\P'l + P'2\ P 

(25) 

P^P'l - P[P'2 



\P!,P'l-P[P'2\ 

By definition, zhx = while zcs can now be expressed as cost i\\ + sinr r being 
the angle between the two axes: 



cos r 



sm r 



icose 

7 



cos^ e + sin^e "^tP 
sm c) Ep'Y 



7^ 



cos2e + sin2e "^tP 
We see that the result depends only on the momentum and mass of the quarkonium 



state (m-T = a/ rn^ + p\)- The angle 5 entering in Eq. 21, equal to r in magnitude 
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defines the positive rotation (respecting the right-hand rule) from one frame to the 
other. Its sign depends, therefore, on the exact conventions used for the orientation 
of the axes y and z of the polarization frames. In the convention where the y axis is 
defined as 

y = ^ ^ 27 

and the z axis is defined by Eq. [25} with the "first" beam oriented as the laboratory 
z axis, the positive rotation is the one bringing the HX axis to coincide with the CS 



axis. We thus write, using Eq. 26 



JHX^-CS — — OCS-s>HX 



arccos . [2t 

uitP 



Equation 21 , containing terms of the form 



sin^ (5hx^cs = sin^ ^cs^hx = ' 



sin25HX-^cs = -sin25cs^HX 



p"^ m\ 
2mpT:pi^E 



(29) 



p2 m\ 



is now explicitly seen as a kinematic-dependent transformation. 

As an example, we show in Fig. [T] how natural i/if) polarizations = +1 and 
— 1 in the CS frame (with A<^ = A^^ = and no intrinsic kinematic dependence) 
translate into different px-dependent polarizations measured in the HX frame in dif- 
ferent rapidity acceptance windows, representative of the acceptance ranges of several 
Tevatron and LHC experiments. Corresponding figures for the T(IS') case can be 
seen in Ref. [21]. The same results, except for a change in the sign of A^<^, the only 
parameter depending on the sign of the rotation angle, are obtained if the roles of 
the two frames are inter-exchanged. 

The change of sign of the rapidity does not change the A^ and A<^ curves. However, 
the sign of A^<^ can change from positive to negative rapidity, depending on the 



convention used for the orientation of the axes. If the axes are defined as in Eqs. 25 



and 27 at both positive and negative rapidity, always taking as "first" beam the one 
positively oriented in the laboratory, A^^, (proportional to sin 25 with < 5 < vr) is 
forced to change sign when the rapidity changes sign. Any measurement integrating 
events over a range in rapidity where the acceptance is symmetrical around zero 
would, therefore, yield A^^ = 0. In order to avoid this cancellation, the axis definitions 



in Eqs. 25 and 27 can be improved by inverting the orientations of y and 5cs for 
negative rapidity (correspondingly restricting the domain of the rotation angle to 
< 5 < 7r/2). 

Having seen how the strong kinematic dependence induced by the choice of the 
observation frame can mimic and/or mask the fundamental ("intrinsic") dependen- 
cies reflecting the production mechanisms, let us now discuss the additional problems 
caused by common experimental limitations. Experiments can only measure the net 
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Figure 7: Kinematic dependence of the ^/ip decay angular distribution seen in the HX 
frame, for natural polarizations = +1 (a-b-c) and = — 1 (d-e-f) in the CS frame. 
The curves correspond to different rapidity intervals; from the solid line: \y\ < 0.6 
(CDF), \y\ < 0.9 (ALICE), \y\ < 1.8 (DO), \y\ < 2.5 (ATLAS and CMS), 2 < |y| < 5 
(LHCb). For simplicity the event populations were generated fiat in rapidity. The 
sign of A^^ depends on the definition of the y axis of the polarization frame, here 
taken as sign(]9L)(-P'i x P'2)/\P'i x P'2\, where are the momenta of the colliding 
particles in the meson's rest frame. 

polarization of the specific cocktail of quarkonium events accepted by the detector, 

trigger and analysis cuts. Moreover, each measured value necessarily implies an inte- 
gration over certain ranges (bins or cells) of the quarkonium momentum components. 
If the polarization depends on the kinematics, this binning effectively biases the 
measured angular parameters, in different ways for experiments having different dif- 
ferential acceptances. In other words, two experiments covering the same kinematic 
interval can measure different average polarizations. This problem can be solved by 
presenting the results in narrow intervals of the probed phase space. Similarly, the- 
oretical calculations aimed at comparisons with experimental data should consider 
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how the momentum distributions are distorted by the acceptances of those exper- 
iments. Ahernatively, the predictions should avoid kinematic integrations or, even 
better, be provided as event-level information to be embedded in the Monte Carlo 
simulations of the experiments. These considerations provide a further motivation 
for reporting measurements and theoretical calculations in frame-independent terms, 
as we will discuss in the next section. 



6 A frame- invariant approach 

The general frame-transformation relations in Eg. [2T] imply the existence of an invari- 
ant quantity, definable in terms of A^, and A^^, in one of the following equivalent 
forms: 

^ (3 + A^) + ci(l-A^) 

C2(3 + A,) + C3(l-A^)- ^^''^ 

An account of the fundamental meaning of the frame-invariance of these quantities 
can be found in Ref. [25]. We will consider here, specifically, the form 

A = ^{-3,0,1} = _ ^ — • I'^lJ 

In the special case when the observed distribution is the superposition of n "elemen- 
tary" distributions of the kind 1 + A^*^ cos^-i?, with event weights f^''\ with respect to 
n different polarization axes, A represents a weighted average of the n polarizations, 
insensitive to the orientations of the corresponding axes: 

A - V A« / V (32) 

i=l ^ + ■^^ I i=l <3 + ^tf 

The determination of an invariant quantity is immune to "extrinsic" kinematic de- 
pendencies induced by the observation perspective and is, therefore, less acceptance- 
dependent than the standard anisotropy parameters A^, A<^ and A^i^. 

This is shown in Fig. |8| where we consider, for illustration, that 60% of the J/?/) 
events have natural polarization A.,? = +1 in the CS frame while the remaining fraction 
has A^ = +1 in the HX frame. Although the polarizations of the two event subsamples 
are intrinsically independent of the production kinematics, in neither frame, CS or 
HX, will measurements performed in different transverse and longitudinal momenta 
windows find identical results. However, in this well as in the simpler case of 

Fig. [7](a-b-c), any arbitrary choice of the experimental observation frame will always 
yield the value A = +1, independently of kinematics. This particular case, where all 
contributing processes are transversely polarized, is formally equivalent to the Lam- 
Tung relation [15], as discussed in Ref. [25]. Analogously, the example represented 
in Fig. [7](d-e-f), or any other case where all polarizations are longitudinal, yields 
A = -1. 
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Figure 8: Kinematic dependence of the i/ip decay angular distribution seen in the 
HX (a-b-c) and CS (d-e-f ) frames, when 60% (40%) of the events have full transverse 
polarization in the CS (HX) frame. The curves represent measurements in different 
acceptance ranges, as detailed in Fig. [7j Corresponding figures for the T(IS') case 
can be seen in Ref. |24j . 



The existence of frame-invariant parameters also provides a useful instrument for 
experimental analyses. Checking, for example, that the same value of an invari- 
ant quantity (Eq. 30) is obtained, within systematic uncertainties, in two distinct 
polarization frames is a non-trivial verification of the absence of unaccounted sys- 
tematic effects. In fact, detector geometry and/or data selection constraints may 
strongly polarize the reconstructed dilepton events. Background processes also affect 
the measured polarization, if not well subtracted. The spurious anisotropies induced 
by detector effects and background do not obey the frame transformation rules char- 
acteristic of a physical J = 1 state. If not well corrected and subtracted, these effects 
will distort the shape of the measured decay distribution differently in different po- 
larization frames. In particular, they will violate the frame-independent relations 
between the angular parameters. Any two physical polarization axes (defined in the 
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rest frame of the meson and belonging to the production plane) may be chosen to 
perform these "sanity tests". The HX and CS frames are ideal choices at high p^, 



where they tend to be orthogonal to each other (in Eq. 29, sin 6 ^ 1 for ^ m 



and/or p^ ^ \Pl\)- At low where the difference between the two frames vanishes, 
any of the two and its exact orthogonal may be used to maximize the significance of 
the test. Given that A is "homogeneous" to the anisotropy parameters, the difference 
A^^-* — A*-^^ between the results obtained in two frames provides a direct evaluation 
of the level of systematic errors not accounted in the analysis. 



7 Effect of intrinsic parton transverse momentum 

In this section we describe how the geometry of the CS frame is related to the kine- 
matics of the production process. It can be recognized from Eq. |26] that the vector 
zcs indicates the direction of the laboratory z axis (that is, the beam line) as seen 
in the quarkonium rest frame. In this frame, any length will be Lorentz contracted 
by a factor I/7 along the quarkonium boost direction, but not along the transverse 
directions. In the quarkonium rest frame (as well as in the laboratory) the direction 
of the beam line coincides with the direction of the relative motion of the colliding 
partons ("parton axis"), when their transverse momenta are neglected (and exactly 
when averaging a large sample of events). This approximation affects the experimen- 
tal determination of an angular distribution naturally of the kind 1 + A^ cos^ 'd* with 
respect to the parton axis z*. In the following considerations we fix a coordinate sys- 
tem having the z axis along the dilepton direction in the laboratory and the xz plane 
coinciding with the production plane. We then define the directions of the beam axis 
and of the parton axis in the laboratory as, respectively, 

B = (sin 9,0, cos 6) , 

. (33) 
B' = (sin 6' cos sin 6' sin cos 6') , 

where 6 and G' are the angles they form with respect to the dilepton direction. The 
presence of the angle $' denotes the fact that, due to the intrinsic transverse momenta 
of the partons, the vector B' does not belong, in general, to the production plane. 
The angle A between the two directions in the laboratory is given by 

cos A = sin B sin 0' cos $' + cos G cos G' . (34) 

When boosted to the dilepton rest frame, the two vectors become 

(sinG,0, ^cosG) 



b' 



sin^ G + cos^ G 
sin G' cos sin G' sin - cos G' 



, (35) 



sin^ G' + ^ cos2 G' 
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and the cosine of the angle between them 



^ cos A — cos O cos 6' 

cosC= — , — , (36) 
v/1 - /32 cos2 e v^l - cos2 6' 

The rotation by ( from the parton axis to the beam hne axis transforms the polar- 
ization parameter X,^ according to the following expressions: 



2 ' - '7 (37) 



These transformations do not represent a simple rotation as Eq. 21 Indeed, the 
magnitude of the polar anisotropy decreases, while no significant azimuthal anisotropy 
arises. In fact, the rotation plane (formed by the parton and beam lines) does not 
coincide with the experimentally defined production plane. The angle between the 
two planes changes from one event to the next, so that the azimuthal anisotropy 
deriving from the tilt between the "natural" polarization axis and the experimental 
axis tends to be smeared out in the integration over all events. Since cos A ~ 1 — 
I sin^ A and (approximately event-by-event, and exactly on average) cos0' ~ cos0, 
from Eq. |36] we obtain 

2 



9 ^> (sin^ A) E , ^ , , ^ 



/32 cos^ m 



Denoting by fci.2, A;i 2t and Ei^2 the total momenta, transverse momenta and energies 
of the two partons in the laboratory, the laboratory angle A satisfies 

sin^A = ^ fciM! (39) 



and, on average. 



where we have defined the average parton squared transverse momentum as {k\) = 
{{k!^) + (4\))/2. 

Considering now the specific case of Drell-Yan production at low p^, we can 
assume an approximate equality between total parton energy and dilepton energy, 
E1 + E2 ~ E, and, moreover, (m^) — m'^ + 2{k'^). Combining Eqs.[37](with AJ = +1) 



38 and 40, we find that the measurement of the polarization of low-px Drell-Yan 



dileptons provides an estimate of the "effective" parton transverse momentum: 
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The average Drell-Yan polarization A^^ = 1.008 ± 0.026 measured by E866 pLSj, in 
proton-copper collisions for (m) ~ 10 GeV/c^ and pt < 4 GeV/c, implies, therefore, 
that 

(kl) < 0.5 GeVVc^ at 68% C.L. and 
J (42) 
(4) < 1.0 GeVVc^ at 95% C.L. . 

Tighter limits could be derived from precise low-mass measurements, given that the 
polarization smearing effects are essentially proportional to m~^. Unfortunately, the 
existing (pion-induced) measurements |26l |27j, though very precise and extending 
down to 4 GeV/c^, present large azimuthal anisotropics of dubious interpretation 
and are scarcely suitable for this purpose. 

We can now estimate the maximum magnitude of the smearing effects that can 
be foreseen for the observable quarkonium polarization when the natural polarization 



axis is the parton axis. Combining again Eqs. 37, [38] and 40, this time with E1 + E2 > 



E, we find that the magnitude of the polarization is reduced by the fraction 



\CS \* 



< 



(3 + AS) (4) 



(43) 



For example, it cannot be excluded, considering the limit in Eq. |42| that a fully trans- 
verse natural polarization of the J/tp along the parton axis is reduced by as much as 
30%, for pt = 2 GeV/c, when observed in the CS frame. This smearing effect should 
be one order of magnitude smaller for J/ip mesons of px = 10 GeV/c. On the other 
hand, given the strong dependence of the effect on the dilepton mass, bottomonium 
polarization measurements are practically insensitive to the parton transverse mo- 
mentum even at low p^- This prediction is consistent with the already quoted E866 
results, showing T(2S' + 35") polarizations in the CS frame always compatible with 
+1, within a ~ 15% uncertainty, in four bins between and 4 GeV/c. 



8 A few concrete examples 

We conclude with some examples of measurements illustrating concepts described in 
the previous sections. 

We have already referred to the E866 measurement of a full transverse T(25' + 3S') 
polarization in the CS frame. The result is represented in Fig. [9]a, as a function of 
Pt- a similarly constant behaviour, consistent with a Drell-Yan-like polarization, has 
been measured by this experiment as a function of xp, in the range [0, 0.5], confirming 
that the adoption of the CS frame is, in this case, an optimal choice. It is true that, 
in special kinematic conditions, the transverse polarization observed in one frame 
could, in reality, be the reflection of a natural longitudinal polarization in another 
frame, as shown in Fig. [7]d. However, a maximal polarization independent of the 
production kinematics in the CS frame must directly reflect the spin configuration of 
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the interacting partons (as is well known to be the case in Drell-Yan production, a 
paradigmatic example of natural transverse polarization in the CS frame). 

To better illustrate the importance of an optimal choice of the reference frame, we 
will now consider what the E866 experiment would have measured, had the analysis 
been made with a different choice. As a reasonable approximation, we assume that the 
azimuthal distribution is exactly isotropic in the CS frame. The polar anisotropy that 
would be observed in the HX frame is shown in Fig. |9]b, where the curve includes an 
extrapolation to higher px assuming that in the CS frame the distribution continues 
to have the shape 1 + cos^ i). A measurement performed in the HX frame would 
show, quite misleadingly, a polarization changing from fully transverse to partially 
longitudinal. The strong signature of "natural" transverse polarization evidenced 
by the data in the CS frame would become unrecognizable, although it could, in 
principle, be reconstructed back if (and only if) the azimuthal anisotropy were also 
measured. 




Figure 9: The E866 measurement of the T{2S + 35*) polarization in the CS (a) and 
a guess of how it would be observed in the HX frame, extrapolating to higher (b). 



Seeing how the curve in Fig. [9]b qualitatively resembles the pattern measured by 
CDF for the J/ip, it is natural to wonder how that measurement, made in the HX 
frame, would look like in the CS frame. Unfortunately, in this case the measure- 
ment itself, a slight longitudinal polarization, does not suggest any educated guess 
on what we could assume for the unmeasured azimuthal anisotropy. For example. 



as shown in Fig. [10} if the distribution in the HX frame were azimuthally isotropic, 
the measured polarization would correspond to a practically undetectable polariza- 
tion in the CS frame (dashed line). However, if we take into account all physically 
possible values of the azimuthal anisotropy, as allowed by the triangular relation rep- 
resented in Fig. |6j we can only derive a broad spectrum of possible CS polarizations, 
approximately included between —0.5 and +1 (shaded band). This example shows 
how a measurement reporting only the polar anisotropy is amenable to several inter- 
pretations in fundamental terms, often corresponding to drastically different physical 
cases. One possible hypothesis would be that all processes are naturally polarized in 
the HX frame and that transverse and longitudinal polarizations are superimposed 
in proportions varying from approximately 2/3 transverse and 1/3 longitudinal at 
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Figure 10: The CDF J ftp polarization measurement in the hehcity frame (data points) 
and the range for the corresponding polarization in the CS, allowing for all possible 
values of the azimuthal anisotropy (shaded band). The dashed line is the CS polar- 
ization for A^-^ = 0. 



Pt = 5 GeV/c (A^ ~ 0) to around 60% /40% at pt = 20 GeV/c (A^ ~ -0.2). In this 
case, no azimuthal anisotropy should be observed in the HX frame. Alternatively, 
we can consider a scenario where the observed slightly longitudinal HX polarization 
is actually the result of a mixture of two processes, both producing J/ip mesons with 
fully transverse polarizations, but one in the HX frame and the other in the CS 
frame. Figure [TTj (left) shows that this scenario is perfectly compatible with the CDF 
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Figure 11: The pt dependence of the angular parameters A^ (left) and A^ (right) 
as would be measured by CDF in the HX frame, according to a scenario where the 
J/V^'s have always full transverse polarization, either in the CS frame or in the HX 
frame, with a suitable px-dependent proportion between the two event samples. The 
data points represent the CDF measurement. 

\^ measurement if the proportion /hx/ (/hx + /cs) between the two sub-processes is 
assumed to vary linearly between 30% at = 5 GeV/c and 15% at pt = 20 GeV/c. 
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The difference with respect to the previous hypothesis is that now we would mea- 
sure a significant azimuthal anisotropy, A<^ ~ 0.3, as shown in Fig. 11 (right). As 
an attempt to reconcile low-px measurements with collider data, Ref. [10] described 
one further conjecture, in which the polarization arises naturally in the CS frame, 
and becomes increasingly transverse with increasing total ^/ip momentum. Again, a 
direct measurement of A<^ (which, in this case, should be zero in the CS frame but 
positive and increasing in the HX frame) would easily clarify the situation. 
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Figure 12: Example of a "gedankenexperiment" where the ^/ip polarization measure- 
ments in the CS and HX frames (empty and filled symbols, respectively) would be 
inconsistent with each other. 

We finish this section by illustrating the application of the frame-independent 
formalism as a tool to estimate residual systematic uncertainties in experimental data 



analyses. Figure 12 shows a putative set oi^/ip polarization measurements performed 
in the CS and HX frames, versus px- While the values seem to change significantly 
from one frame to the other, the two A;^ patterns are very similar. This observation 
should alert to a possible experimental artifact in the data analysis. We can evaluate 
the significance of the apparent contradiction by calculating the frame-invariant A 
variable in each of the two frames. For the case illustrated in Fig. [121 averaging the 



four represented bins, we see that A in the HX frame is larger than in the CS 
frame by 0.5 (a rather large value, considering that the decay parameters are bound 
between —1 and -|-1). In other words, an experiment obtaining such measurements 
would learn from this simple exercise that its determination of the decay parameters 
must be biased by systematic errors of roughly this magnitude. Given the puzzles and 
contradictions existing in the published experimental results, as recalled in Section [1} 
the use of a frame-invariant approach to perform self-consistency checks, which can 
expose unaccounted systematic effects due to detector limitations and analysis biases, 
constitutes a non-trivial complementary aspect of the new approach proposed in this 
paper for quarkonium polarization measurements. 
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9 Summary and conclusions 



Motivated by several puzzles affecting the existing measurements of quarkonium po- 
larization, we present in this paper a set of proposals which should improve the 
experimental determination of the ^/ip and T polarizations. They are summarized 
in the next paragraphs. 

Measurements and calculations of vector quarkonium polarization should provide 
results for the full dilepton decay angular distribution (a three-parameter function) 
and not only for the polar anisotropy parameter. Only in this way can the measure- 
ments and calculations represent unambiguous determinations of the average angular 
momentum composition of the produced quarkonium state in terms of the three base 
eigenstates, with Jz = +1, 0, —1. 

It is advisable to perform the experimental analyses in at least two different polar- 
ization frames. In fact, the self-evidence of certain signature polarization cases (e.g. 
a full polarization with respect to a specific axis) can be spoiled by an unfortunate 
choice of the reference frame, which can lead to artificial ("extrinsic") dependencies 
of the results on the kinematics and on the experimental acceptance. 

The measured dependence of the polarization on the production kinematics is 
necessarily influenced by the differential experimental acceptance, i.e. by the kine- 
matic distribution of the population of the accepted events. This problem, which is 
not solved by acceptance corrections, can be minimized by providing the results in 
narrow cells in quarkonium rapidity and transverse momentum. Theoretical calcula- 
tions should be provided as event-level inputs to Monte Carlo generators which can 
be tailored to the specific performance capabilities of each experiment. 

The decay angular distribution can be characterized by a frame-independent quan- 
tity, such as A, calculable in terms of the polar and azimuthal anisotropy parameters. 
The existence of such frame-invariant quantities can be used during the data analy- 
sis phase to perform self-consistency checks that can expose previously unaccounted 
biases, caused, for instance, by the detector limitations or by the event selection 
criteria. 

Besides providing a much needed control over systematic experimental biases, the 
variable A also provides relevant physical information: it characterizes the shape of 
the angular distribution, reflecting "intrinsic" spin- alignment properties of the decay- 
ing state, irrespectively of the specific geometrical framework chosen by the observer. 
For instance, we obtain A = +1 for the shapes shown in Fig. [s] (a) and (c), and 
A = — 1 for the shapes shown in the panels (b) and (d). A very important advantage 
of re-expressing the frame-dependent polar and azimuthal anisotropics in terms of a 
frame-invariant quantity is the exact cancellation of extrinsic dependencies on kine- 
matics and acceptances, enabling more robust comparisons with other experiments 
and with theory. The calculation of A requires, anyhow, the determination of the full 
decay distribution in a chosen reference frame and, obviously, does not replace this 
standard procedure. Moreover, the three frame-dependent parameters (A^, A<^ and 
A^ip) can provide information on the direction of the spin- alignment of the decaying 
particle (when this direction is univocally defined) and, therefore, on the topological 
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properties of the dominant production mechanism. For instance, the measurement 
of a full transverse polarization in the CS frame represents a direct observation of 
the spin alignments of the interacting partons, as we know from Drell-Yan produc- 
tion. On the other hand, in the presence of a superposition of production processes 
with polarizations along different axes, measuring the frame-dependent anisotropies 
will not provide, in general, much information on the polarizations involved or on 
their natural alignment directions, while the value of A will immediately tell us if the 
processes involved have a predominantly transverse or longitudinal nature. 

Stripped-down analyses which only measure the polar anisotropy in a single ref- 
erence frame, as often done in past experiments, give more information about the 
frame selected by the analyst ("is the adopted quantization direction an optimal 
choice?") than about the physical properties of the produced quarkonium ("along 
which direction is the spin aligned, on average?"). For example, a natural longitu- 
dinal polarization will give any desired value, from —1 to +1, if observed from 
a suitably chosen reference frame. Lack of statistics is not a reason to "reduce the 
number of free parameters" if the resulting measurements become ambiguous. The 
forthcoming measurements of quarkonium polarization in proton-proton collisions at 
the LHC have the potential of providing a very important step forward in our under- 
standing of quarkonium production, if the experiments adopt a more robust analysis 
framework, incorporating the ideas presented in this paper. 
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